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The differential cross section for the process 77 titt^ has been measured in the kinematic 
range 0.84 GeV < W < 4.0 GeV, |cosS*| < 0.8, where W and 0* are the energy and n° (or 
77) scattering angle, respectively, in the 77 center-of-mass system. The results are based on a 
223 fb~^ data sample collected with the Belle detector at the KEKB e~^e~ collider. Clear peaks 
due to the ao(980) and a2(1320) are visible. The differential cross sections are fitted in the energy 
region 0.9 GeV < W < 1.46 GeV to obtain the parameters of the ao(980). Its mass, width and 
rj^B{riTT°) are measured to be 982.3 to ? -1.7 MeV/c^ 75.6 ± 1.6 tlH MeV and 128 tl tlf eV, 
respectively. The energy and angular dependences above 3.1 GeV are compared with those measured 
in the tt'^tt" channel. The integrated cross section over |cos6'*| < 0.8 has a W~" dependence 
with n = 10.5 ± 1.2 ± 0.5, which is slightly larger than that for tt^tt". The differential cross 
sections show a sin"'* ^* dependence similar to 77 Tr^yr". The measured cross section ratio, 
a{rin°)/a{n°TT°) = 0.48 ± 0.05 ± 0.04, is consistent with a QCD-based prediction. 

PACS numbers: 13.60.Le, 13.66.Bc, 14.40. Cs,14.40.Gx 



I. INTRODUCTION 



Measurements of exclusive hadronic final states in two-photon collisions provide valuable 
information concerning physics of light and heavy-quark resonances, perturbative and non- 
perturbative QCD and hadron-production mechanisms. So far, we have measured the pro- 
duction cross sections for charged-pion pairs charged and neutral- kaon pairs 
and proton-antiproton pairs [6|]. We have also analyzed D-meson-pair production and ob- 
serve a new charmonium state identified as the Xc2(2P) 0|- In addition, we have measured 
the 7r°7r° final state [1, The statistics of these measurements is two to three orders of 
magnitude higher than pre-B-factory measurements opening a new era in studies of 
two-photon physics. 

In this paper, we report measurements of the differential cross sections, da/d\ cos 6**1, for 
the process 77 — > t^tt" in a wide two-photon center-of-mass (cm.) energy (W) range from 
0.84 GeV to 4.0 GeV and in the cm. angular range, | cos6'*| < 0.8. We use only the 77 77 
and 7r° — *■ 77 decay modes in this analysis. The decay mode 77 7r"'"7r~7r° is not used 
because of a much lower product of efficiency and branching fraction. 

Previously, it was reported that this reaction is dominated by resonance production 
We can restrict the J^'" quantum numbers of the meson produced by two photons to 
be l~(even)"'""'", that is, those of aj^evon mesons. A long-standing puzzle in QCD is the 
existence and structure of low mass scalar mesons. In the / = sector, we recently observed 
a peak for the /o(980) in both the 77 tt+tt" and 77 7r°7r° channels The two- 

photon width of the /o(980) is measured to be 200 — 300 eV, supporting its nature (l^ . 
Our analysis also suggests the existence of another /o meson in the 1.2-1.5 GeV region that 
couples with two photons @. In the 1=1 sector, the an (980) and a2(1320) have been 
observed previously with a rather low statistical significance [HI, The parameters for 
the ao(980), in particular its two-photon width are of great interest because of its connection 
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to the nature of low mass scalar mesons. Moreover, other scalar or tensor aj mesons can be 
searched for in the higher mass region above the prominent peak from the 02(1320). 

If there were an / = 1 "hidden-charm" (that is, charmonium-like) meson, it could be a very 
strong candidate for an exotic state, because charmonia have / = and isospin is conserved 
in their hadronic deca^. However, recently, some new particles that are not pure / = 0, 
such as the X(3872) [iJI and Z(4430) [Hi, have been reported. 

At higher energies {W > 2.4 GeV), we can invoke a quark model. In leading-order calcula- 
tions, the ratio of the rjir^ cross section to that of tt^tt" is predicted within uncertainties due 
to the different form factors for the 7r° and rj. Analyses of energy and angular distributions 
of these cross sections are essential to determine properties of the observed resonances and 
to test the validity of QCD based models [l^.[r^. 

This paper is organized as follows. In Sec. [ill the experimental apparatus used and the 
event selection are described. Section IIIII explains background subtraction and derivation 
of the differential cross sections. In Sec. IIVI the resonance parameters of the ao(980) are 
derived by parameterizing partial wave amplitudes with resonances and smooth non-resonant 
background amplitudes and fitting differential cross sections. Section |V] describes analyses 
at higher energy. The topics included there are the angular dependence of differential cross 
sections, the W dependence of the total cross section, and the ratio of cross sections for rjn^ 
to 7r°7r° production. Finally, Section IVll summarizes the results and presents the conclusion 
of this paper. 

II. EXPERIMENTAL APPARATUS AND EVENT SELECTION 

Events consisting only of neutral final states are extracted from the data collected in 
the Belle experiment. In this section, the Belle detector and event selection procedure are 
described. 



A. Experimental apparatus 

A comprehensive description of the Belle detector is given elsewhere • We mention here 
only those detector components that are essential for the present measurement. Charged 
tracks are reconstructed from hit information in the silicon vertex detector and the central 
drift chamber (CDC) located in a uniform 1.5 T solenoidal magnetic field. The detector 
solenoid is oriented along the z axis, which points in the direction opposite to that of the 
positron beam. Photon detection and energy measurements are performed with a CsI(Tl) 
electromagnetic calorimeter (ECL). 

For this all-neutral final state, we require that there be no reconstructed tracks coming 
from the vicinity of the nominal collision point. Therefore, the CDC is used for vetoing 
events with charged track(s). The photons from decays of the neutral pion and the rj meson 
are detected and their momentum vectors are measured by the ECL. The ECL is also used to 
trigger signal events. Two kinds of the ECL trigger are used to select events of interest: the 
ECL total energy deposit in the triggerable acceptance region (see the next subsection) is 
greater than 1.15 GeV (the "HiE" trigger), or the number of ECL clusters counted according 
to the energy threshold at 110 MeV for segments of the ECL is four or larger (the "Clst4" 
trigger). The above energy thresholds are determined by studying the correlations between 
the two triggers in the experimental data. No software filtering is applied for triggering 
events by either or both of the two ECL triggers. 
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TABLE I: Data sample: luminosities and energies 



e'^e cm. energy 


Luminosity 


Runs 


(GeV) 


(fb-^) 




10.58 


179 


T(4S) 


10.52 


19 


continuum 


10.36 


2.9 


r{3S) 


10.30 


0.3 


continuum 


10.86 


21.7 


r(5S) 


Total 


223 





B. Experimental data and data filtering 

We use a 223 fb^^ data sample from the Belle experiment at the KEKB asymmetric-energy 
e+e~ collider [l^. The data were recorded at several e+e~ cm. energies summarized 
in Table HI The difference of the luminosity functions (two-photon flux per e+e^-beam 
luminosity) in the measured W regions due to the difference of the beam energies is small 
(maximum ± 4%). We combine the results from the different beam energies. The effect on 
the cross section is less than 0.5%. 



The analysis is carried out in the "zero-tag" mode, where neither the recoil electron nor 
positron are detected. We restrict the virtuality of the incident photons to be small by 
imposing a strict requirement on the transverse-momentum balance with respect to the 
beam axis for the final-state hadronic system. 

The filtered data sample ("Neutral Skim") used for this analysis is the same as the one 
used for tt^tt^ studies The important criteria in this filtering are: no good tracks; two 

or more photons or one or more neutral pions that satisfy a specified energy or transverse- 
momentum criterion. Performance of the ECL triggers is studied in detail using the 7r°7r° 
events 

C. Event selection 

From the Neutral Skim event sample, we select 77 rjir^ with the following conditions: 

(1) the total energy deposit in ECL is smaller than 5.7 GeV; 

(2) there are exactly four photons in the ECL each having energy larger than 100 MeV; 

(3) the ECL energy sum within the triggerable region is larger than 1.25 GeV, or, all the 
four photons are within the triggerable region, i.e. in the polar-angle range, —0.6255 < 
cosO < -1-0.9563, in the laboratory frame; 

(4) a combination of two photons is reconstructed as a neutral pion that satisfies the 
following conditions on invariant mass, |A'/(77) - 0.1350 GeV/c^| < 0.0200 GeV/c^, 
transverse momentum pt(7r°) > 0.15 GeV/c and goodness of the mass-constrained fit 

(5) the combination of the remaining two photons has an invariant mass consistent with 
r] 77, 0.51 GcV/c^ < M(77) < 0.57 GeV/c^. There are three combinations of 
photon pairs that can be constructed from the four photons, and all the combinations 
are tried and any of them satisfying the above criteria are retained. We scale the 
energy of the two photons in (5) with a factor that is the ratio of the nominal 77 mass 
to the reconstructed mass. This is equivalent to an approximate 1-C mass constraint 
fit where the relative energy resolution (AE/E) is independent of E and the resolution 
in the angle measurement is much better than that of the energy. The present case 
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PIG. 1: Two-dimensional W and |cos0*| distribution of the r/n^ candidates in data. The same 
distribution is viewed from two different directions. 

is close to this. After scaling the 77's four-momentum, we calculate the invariant mass 
(W) and the transverse momentum |) in the e~^e~ cm. frame for the rjir'^ system. 
(6) The transverse momentum is is required to be less than 0.05 GcV/c. 

We define the cm. scattering angle, 0*, as the scattering angle of the 7r° (or equivalently, 
of the rj ) in the 77 cm. frame, for each event. We use an approximation that the e+e~ 
axis is the reference for this polar angle as (since we do not know the exact 77 axis in the 
zero-tag condition). 

Signal and background events for e"'"e~ e^e~r/Tr'^ are generated using the TREPS 
code [iOl. All Monte Carlo (MC) events are put through the trigger and detector simulators 
and the event selection program. We find that up to 2% of events in the region below 
W < 1.05 GeV have two entries per event because of the multiple combinations satisfying 
criteria (4), (5) and (6). The two entries per event have similar W and |, but different 
I cos 6** I values. The fraction with double entries is small and is in principle compensated 
by the normalization using the efficiency determined by the MC sample. We find a similar 
fraction of multiple entries in the signal-MC data. 

A total of 2.82 x 10^ events are selected from 3.53 x 10* events of the Neutral Skim 
sample. The lego plots of two-dimensional distributions of the selected events (after requiring 
lEpt I < 50 MeV/c) are shown in Fig. [TJ The projected W distribution integrated over 
I 008 6**1 < 0.8 is shown in Fig.O We find at least three resonant structures: near 0.98 GeV 
(ao(980)), 1.32 GeV (a2(1320)) and 1.7 GeV (probably the a2(1700)). 

III. DERIVING DIFFERENTIAL CROSS SECTIONS 

In this section, we present the procedure to derive differential cross sections. First, the 
nature and origin of backgrounds and the method for their subtraction are described. Un- 
folding is then appHed, efficiencies are determined, differential cross sections are derived and 
their systematic errors are estimated. 
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FIG. 2: W distribution for rjiv'^ candidate events integrated over |cos^*| < 0.8. The solid curve is 
the pt-unbalanced background that are experimentally determined. The dashed curve includes an 
additional contribution from the rj' TyTr'^yr" background. 



A. Background subtraction 

In the entire energy range the background is primarily from photons originated from spent 
electrons as identified by unbalanced pt. At low energy there is additional background from 
the rj' decaying into rjir^TT^. 



1. Background from r( titt^tv'^ 

A scatter plot of pt balance vs. invariant mass for the candidate events, Fig. [Slja), shows 
a concentration of events in the pt-unbalanced region in the vicinity of 0.82 GeV. This 
energy is close to the mass difference between the rj' and 7r°, and this structure is due to 
the background from rj' r/n'-'Tr'^ where two photons from one 7r° are undetected. Since 
the background is larger than the signal in the r]' mass region, we cannot measure the cross 
section below W < 0.84 GeV. Kinematically, for rj' decays the invariant mass of the detected 
7771° system cannot be greater than 0.823 GcV/c^. However, there is a rather long tail on 
the higher side up to ~ 0.88 GeV. This background is due to a fake pion reconstructed 
from a photon with another photon from a 7r° or a noise cluster. Figure (Sljb) shows that t] 
mesons are relatively cleanly reconstructed. 

We subtract the background from primary 77's in two-photon collisions using MC. The 
normalization of the 77' — *■ •qn^n^ background is determined using the experimental data in 
the control region, 0.80 GeV <W< 0.84 GeV integrated over all angles and the transverse- 
momentum range, < 0.15 GeV/c (Fig. |4lja)). We assume that the signal yield is 
negligibly small in this W range. The pt distribution is decomposed into the rj' background 
peaking near 0.05 GcV/c and the other pt-unbalanced component by performing a fit, where 
we use the signal and background functional shapes described in the next subsection for the 
77' and the latter background, respectively. The shape parameters of the rj' component are 
fixed to those obtained from the fit to the corresponding MC data. 

Prom the product of the two-photon decay width and the branching fraction of the 77', 
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FIG. 3: (a) Two-dimensional {W, \ '^Pt\) distribution of the rj-K^' candidate events in the relatively 
low-W region. (b)Two-dimensional distribution of the invariant mass of the two detected photons 
near the rj mass and the Pt-balance of events calculated from the mass of four photons (before 
applying the mass constraints) for events with M{r]'iT°) < 1.0 GeV/c^. 



^■yy{v')^{v' f?7r°7r°), we can estimate the absolute size of the background yield from this 
source. The ratio of the observed background to the MC expectation is 0.91 ± 0.02{stat.) ± 
0.09{syst.) ± 0.06(r.y^S), where the first two errors are experimental, and the last error is 
from the uncertainty of the known 77' properties. This factor is consistent with unity. 

Using the normalization thus determined and the MC events, the background yields from 
this source in each angular bin in the range 0.84 GeV < W < 0.90 GcV are determined. 
This background component is incorporated in the fit described in the next subsection with 
the yield and shape fixed. We neglect the 77' background in the W region above 0.90 GeV. 



2. Subtraction of backgrounds 

In the W region below 2.0 GeV, the pt-unbalanced component is non-negligible and is 
subtracted by fitting the pt distributions. The fitting function is a sum of the signal and 
background components. The 77' background estimated in the previous subsection is incor- 
porated in the fit with a fixed shape and size for W < 0.90 GeV (Figs. |4][b-d)). 

The signal component follows an empirical parameterization from signal MC: 

Ax 

{x is \ '}2pI\, A, B and C are fitting parameters), where the distribution has a linear shape 
near a; = and decreases as ^ x^^'^ at large x. Here B is an important parameter that 
determines the peak position; the peak is at a; = B^ /l.l. The background is parameterized 
by a linear function vanishing at x = in a; < 50 MeV/c and a second-order polynomial for 
a; > 50 MeV/c connected smoothly (up to the first derivative) at a; = 50 MeV/c. Fits are 
applied for < 200 MeV/c in each bin of 1^ < 2.0 GeV and | cos 6**1 < 0.8 with the bin 

widths of 0.04 GeV and 0.1 for the two directions, respectively. 

The background yields found from the fits are fitted to a smooth two-dimensional function 
of (W ^ I cos 0*1), in order to minimize the statistical fluctuations from the MC simulation. 
In Fig. O the curves corresponding to the background thus determined are shown, as well 
as the fixed background from 77' decays. 

The backgrounds are subtracted from the experimental yield distribution. The 77' back- 
ground estimated in the previous subsection is also subtracted aXW < 0.90 GeV. The error 
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FIG. 4: (a): The pt distribution of the rjTv candidates at W = 0.80-0.84 GeV integrated over 
I cos 5* I < 0.8. The peak near 0.05 GeV/c is attributed to the background from rj' — > TyTr'^Tr" and 
the normalization is determined by a fit to MC events of this process. The dashed curve shows 
the experimentally determined background from the pt-unbalanced background component and the 
solid curve is the sum of the two background components, (b), (c) and (d): The pt distribution of 
the 7771" candidate events in the shown kinematical regions. Estimated background yields of the two 
components are shown by the curves (explained in the caption for (a)). The excess over the solid 
curve near j X] P* I ~ corresponds to the signal process. 



arising from this background subtraction is taken into account as a systematic error (see the 
subsection E). 

To confirm the vaHdity of the background subtraction method, we examine the W depen- 
dence of the yield ratio R defined as: 

^ y(0.15 GeV/c < I < 0-20 GeV/c) 

niEKl <0.05GeV/c) ' ^ ' 

where Y is the yield in the specified region. We observe that structures have similar features 
as those seen in the estimated pt-unbalanced background (Fig. [5]). There is a peaking 
structure in the background around W = 1.5 GeV, and a small enhancement just above 
2.0 GeV. No significant structures are seen from 2.5 GeV up to ~ 3.3 GeV. 

We observe a significant increase of R above 3.3 GeV. Such an increase is not reproduced 
by the MC simulation, where only a slow increase of R {R = 0.15 at 2.0 GeV and R = 0.25 at 
4.0 GeV) is expected. The excess of R, AR, is the contribution from multibody background 
processes. We apply a background subtraction with the background contamination level 
estimated to be AR/6. This factor is obtained by assuming a quasi-Hnear | Pt I dependence 
of the background and extracting its leakage into the signal region (| X]Pt I < O-O^ GeV/c), 
which is approximately 1/6 of the yield in the 0.15 GeV/c < | X)Pt I < 0-20 GeV /c region. 
We take a half of the correction (i.e. Ai?/12) as the systematic error from this source if this 
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FIG. 5: Energy dependence of the ratio of the yield in the pt-unbalanced region to that in the 
balanced region (see text for the exact definition), which indicates the level of background contam- 
ination. 

systematic is larger than the 3% uncertainty nominally applied for the whole region above 
2.0 GeV (see the section on systematic errors). The actual sizes of the correction used (i.e. 
AR/6) in W bins are 3.8% (3.35 GeV), 8.5% (3.65 GeV), and 13.0% (3.95 GeV), respectively. 
These corrections are much smaller than the statistical errors in these W bins. 

We find that the background in the rj mass sideband is negligibly small after the subtrac- 
tion of these backgrounds; there are no pt-balanced backgrounds in the rj mass sideband of 
the pt-balance distribution shown in Fig. [3ljb). Therefore, we do not perform background 
subtraction for the non-77 component. 



We unfold the experimental yields in the W distributions to correct for the finite invariant- 
mass resolution in the measurement. The smearing function is an asymmetric Gaussian 
function, which is determined by the signal MC with a further empirical correction. The 
standard deviations of the function are assumed to follow 1.3 x (1.4 — 0.3/W^)% {W is in 
GeV, and the resolution varies by 1.3%-1.8% depending on W) on the lower side of the peak 
and 0.77 x (1.4 - 0.3/^^^)% (varies by 0.8%-l.l%) on the higher side. The resolution is 
slightly better than that in the 7r°7r° case in the low W region, because of the large opening 
angle of the two photons from 77 decay. 

The unfolding procedure is applied for 0.9 GeV < < 1.6 GeV with a bin width of 
0.02 GeV, and for 1.6 GeV <W < 2A GeV with a bin width of 0.04 GeV. The unfolding is 
done independently in each angular bin, whose width is A| 008 6**1 = 0.05 for W < 1.6 GeV 
and A| cos 0*1 = 0.1 for W > 1.6 GeV. Figure [6] shows the yield distributions before and 
after the unfolding in the smallest | cos0*| bins. At higher energies, no unfolding is applied 
since the experimental yield is still insufHcient. 



B. Unfolding the W distributions 
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FIG. 6: The yield distributions before (triangles) and after (diamonds) the unfolding in the smallest 
I cos 61*1 bins, (a) | cos 61*1 < 0.05 and (b) |cose*| < 0.1. 



C. Determination of the efficiency 



The signal MC events for e+e^ e^e^rjn^ are generated using the TREPS code [131 for 
the efficiency calculation at 36 fixed W points between 0.75 and 4.2 GeV and isotropically in 
I cos 9* I . The angular distribution at the generator level does not play a role in the efficiency 
determination, because we calculate the efficiencies separately in each | cos 6**1 bin with a 
width of 0.05. The Q,^„ax parameter that gives a maximum virtuality of the incident photons 
is set to 1.0 GeV^, and the form factor for the cross sections for the virtual photon collisions, 
a.y~^{0,Q^) = a~^~f{0,0)/{l + Q^/W^)^ is used. This form factor does not play any essential 
role in the present analysis, since our stringent pt-balance cut (|X]Ptl < O-O^ GeV/c) 
requires Q^/W^ for the selected events to be much smaller than 1. A sample of 400,000 
events is generated at each W point and is subjected to the detector and trigger simulations. 
The obtained efficiencies are fitted to a two-dimensional function of {W, | cos 6**1) with an 
empirical functional form. 

We embed background hit patterns from random trigger data into MC events. We find 
that different samples of background hits give small variations in the selection efficiency 
determination. A W^-dependent error in the efficiency, 2-4%, arises from the uncertainty in 
this effect. Figure [7] shows the two-dimensional dependence of the efficiency on ( W, \ cos 6**1) 
after the fit for smoothing. 



D. Derivation of differential cross sections 



The differential cross section for each (W, \ cos9*\) point is given by: 

da AY -AB 

d|cos6i*| " AW A\ cos 0*\ J CdtL^^{W)r] ' 

where AY and AB are the signal yield and the estimated pt-unbalanced background in the 
bin, AW and A| cos 6**1 are the bin widths, / Cdt and L^^^jW) are the integrated luminosity 
and two-photon luminosity function calculated by TREPS [20|, respectively, and rj is the 
efficiency including the correction described in the previous section. The bin sizes for W 
and A I cos9*\ are summarized in Table Ull 
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FIG. 7: Two-dimensional dependence of the eflEiciency on {W , | cos 6^*1). 



TABLE II: Bin sizes 



W range 


AW 


A| cose* 


(GeV) 


(GeV) 




0.84 - 1.6 


0.02 


0.05 


1.6 - 2.4 


0.04 


0.10 


2.4 - 4.0 


0.10 


0.10 



Figure [8] shows the angular dependence of the differential cross sections for some selected 
W regions. Figures [9] show the cross section integrated over |cos0*| < 0.8 on logarithmic 
and linear scales for partial W regions. The data points are in good agreement with those 
of Crystal Ball 

E. Systematic errors 

Various sources of systematic uncertainties assigned for the rjn^ signal yield, efficiency 
and the cross section evaluation are described in detail below and summarized in Table IIIII 

Trigger efficiency: The systematic error due to the Clst4 trigger (the mnemonics of the 
ECL triggers were described in Sec. II. A) is taken to be 2/3 of the difference in the efficiencies 
when the thresholds for the Clst4 trigger are varied from 110 MeV to 100 MeV in the trigger 
simulator, and no error is assigned for W > 2.5 GeV where the HiE trigger plays a dominant 
role. In addition, we take the uncertainty in the efficiency of the HiE trigger to be 4% in the 
whole W region. The systematic errors from the two triggers are combined in quadrature. 
The former component is approximated by an angular-independent function of the cm. 
energy, 27% x Q.^mw-a.sb [GcV])^ rjj^jg exceeds 5% for W < 1.2 GeV. 

The reconstruction efficiency: we assign 6% for the reconstruction of a tt'^ and an rj. 

Pt-balance cut: 3% is assigned. Thept-balance distribution for the signal is well reproduced 
by MC so that the efficiency is correct to within this error. 
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FIG. 8: Angular dependence of the differential cross sections for eight selected W bins indicated. 
The bin sizes are summarized in Table Ull 



Background subtraction: 20% of the size of the subtracted component is assigned to this 
source for the range 0.9 GeV < W < 2.0 GeV. In the W region where the background 
subtraction is not appHed {W > 2.0 GeV), we assign a systematic error of 3%, which 
is a conservative upper Hmit on the background contamination from an investigation of 
the experimental distributions. Above 3.5 GeV, the error originating from the background 
subtraction (Ai?/12) is larger than 3%, and we replace the error by the latter value. In 
the region 0.84 GcV <W< 0.90 GeV where the 77' background is subtracted, the assigned 
error is a quadratic sum of 7% of the subtracted 77' background and 20% of the subtracted 
P4-unbalanced background. 

Luminosity function: We assign 4% (5%) for W < (>) 3.0 GeV. 

Beam background effect for event selection: We assign a 2% - 6% error depending on 
W for uncertainties of the inefficiency in event selection due to beam-background photons. 
The uncertainty is estimated from the variation of efficiencies among different experimental 
periods or background conditions. We adopt the averaged efficiency from the different 
background files, and the uncertainty in the average is assigned as the error. 

Unfolding: Uncertainties from the unfolding procedure, using the single value decompo- 
sition approach in Ref. [131, are estimated by varying the effective-rank parameter of the 
decomposition within reasonable bounds. 

Other efficiency errors: An error of 4% is assigned for uncertainties in the efficiency 
determination based on MC including the smoothing procedure. 

The total systematic error is obtained by adding all the sources in quadrature and is 10- 
12% for the intermediate and high W regions. It becomes much larger for < 1.06 GeV. 



IV. STUDY OF RESONANCES 



In this section, we extract the resonance parameters of the ao(980) and a possible reso- 
nance ao(1450), as well as check the consistency of the 02(1320) parameters. We also study 
whether or not the 02(1700) is produced in this reaction. 
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FIG. 9: The cross section integrated over |cosS*| < 0.8 on a logarithmic (a) and linear (b) scale 
compared with the Crystal Ball measurement (|cos^*| < 0.9) The corrections for different 

j cos9*\ coverage are not made. The dashed curve shows the size of the systematic error. 



TABLE III: Systematic errors for the differential cross sections. Ranges of errors are shown when 
they depend on W. 



Source Error (%) 



Trigger efficiency 


4 


- 30 


rj and 7r° reconstruction efficiency 




6 


pt-balance cut 




3 


Background subtraction 


3 


- 35 


Luminosity function 


4 


- 5 


Overlapping hits from beam background 


2 


- 6 


Unfolding procedure 





- 4 


Other efficiency errors 




4 



Overall 10 - 12 (for W > 1.06 GeV) 
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FIG. 10: Resulting spectra of 5^ Dq and DI The error bars shown are diagonal statistical errors. 



Formalism 



The formalism is exactly the same as that for tttt P, 0, H, [§| and the analysis is quite 
similar. In the energy region W < 2.0 GeV, J > 2 partial waves (the next is J = 4) may 
be neglected so that only S and D waves are to be considered. The differential cross section 
can be expressed as: 



da , 



(4) 



where S represents the S wave, Dq {D2) denotes the helicity (2) components of the D wave, 
respectively, and Y}" are the spherical harmonics. Since the Idy l's are not independent 
of each other, partial waves cannot be separated using measurements of differential cross 
sections alone. To overcome this problem, we write Eq. ^ as: 



da 



A'n:d\ cost 



,(77 - = \Y,r + Dl + Dl |r,Y 



(5) 



The amplitudes 5^, and correspond to the cases where interference terms are ne- 
glected; they can be expressed in terms of S*, Dq and D2 as follows Q: 



^0 



ISl"" + V5Rc{S*Do) , 
\Do\^ + ■^Rc{S*Do) , 



(6) 



Since squares of spherical harmonics are independent of each other, we can fit the differential 
cross section at each W to obtain S'^, Dq, and D^. The unfolded differential cross sections 
are fitted by taking into account statistical errors only, which will not be independent at 
each W because of the unfolding procedure. However, we treat them as independent in the 
fit. The resulting S*^, Dq and spectra for < 2 GeV are shown in Figs.fTOl 



B. Fitting partial wave amplitudes 



Although the derived amplitudes 5*^, Dq and are functions of partial waves (Eq. ([6])), 
they do give some indication of their behavior. Notably, the Do wave appears to be small 
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and the D2 wave is dominated by the a2(1320) with a hint of the 02(1700). The peak in 
around W = 0.98 GeV (Fig. [T0| is clearly due to the ao(980) resonance and a shoulder 
above the ao(980) peak may be due to the ao(1450). 

In this section, we derive information about resonances by parameterizing partial wave 
amplitudes and fitting differential cross sections. Note that we do not fit the obtained S*^, 
Dq and D2 spectra but fit the differential cross sections directly. Once the functional forms 
of amplitudes are assumed, we can use Eq. ^ to fit differential cross sections. We can then 
neglect the correlations between S*^, Dq and Z?!. The S*^, Dq and -D^ spectra are used to 
make an initial determination of which resonances are important and to check fit quality. 

The energy range of this measurement can be naturally divided into a low energy region 
{W < 1.5 GeV), where the ao(980) and 02(1320) are important, and a higher energy region 
{W > 1.5 GeV), where the parameters of the 02(1700) will be of interest. In both regions, 
we neglect J > 2 waves. 

First we extract parameters of the O2(1700) in the high energy range. A hint of the 
02(1700) is visible in Fig. [9l Various fits have been performed: a fit in the energy region 
1.5 GeV < W < 2.0 GeV with the parameters of the 02 (1700) all fioated or the mass and 
width fixed to the PDG values; a fit in the energy region 0.9 GeV <W< 2.0 GeV by fixing 
the lower energy parameters to the ones determined below and either fioating or partially 
fixing the parameters of the O2(1700). The resulting parameters of the O2(1700) vary much. 
In addition, the fit quality is poor. Thus, we cannot report a definite conclusion on the 
02(1700). 

In the low energy region, W < 1.5 GeV, we can safely neglect J > 2 waves. The Dq 
contribution is small while l)| is seen to be dominated by the O2(1320) resonance. We 
assume that the 02 (1320) contributes to the D2 wave only, since is small. The shoulder 
in the 5^ spectrum above the oo(980) peak may be due to the oo(1450). However, in this 
fit we introduce a new resonance oo(F) instead, since its parameters are found to be quite 
different from those of the oo(1450). The goal of analysis is to obtain parameters of the 
oo(980), 00(F) and to check the consistency of the O2(1320) parameters that have been 
measured well in the past. 

1. Parameterization of amplitudes 

We parameterize S, Dq and D2 waves as follows. 

S = A,„(98o)e'*- + A,„(y)e'^- + Bs , 

D2 = Aa.(1320)e^^^^ +Sd2 , (7) 

where ^ao(980)) ^ao(Y) and ^02(1320) a-re the amplitudes of the Oo(980), oo(F) and 02(1320), 
respectively; I3s, Bjjo and B]j2 are non-resonant (called hereafter "background") amplitudes 
for S, Do and D2 waves; and 4'so, (j)si, and 0^2 are the phases of resonances relative to 
background amplitudes. We also study the case with no 00(F) and the case with the mass 
of the oo(r) fixed to that of the ao(1450). 

The background amplitudes are parameterized as follows. 

Bs = asW'^ + bsW' + Cs + i{a'^W'^ + b'^W' + c'^), 
Bdo = aoW'^ + boW' + co + i(a'oW'^ + b'oW' + c'q), 
Bd2 = a2W'^ + b2W' + C2+i{a'2W'^ + b'2W' + c'2). (8) 

Here W = W — Wth where Wth is the threshold energy. We assume background amplitudes 
to be quadratic in W for the both real and imaginary parts of all waves. In this way, 
symmetries among amplitudes are kept. The arbitrary phases are fixed by choosing ^^o = 
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TABLE IV: Parameters of the a2(1320) 



Parameter 


Value 


Unit 


Mass 


1318.3 ±0.6 


MeV/c^ 


Ptot 


107 ±5 


MeV 


B{a2 pn) 


70.1 ±2.7 


% 


B{a2 777r) 


14.5 ± 1.2 


% 


B{a2 LUTvn) 


10.6 ±3.2 


% 


B{a2 KK) 


4.9 ±0.8 


% 


B{a2 77) 


(9.4 ±0.7) X 10" 


-6 



4>d2 = 0. We constrain all the background amplitudes to be zero at the threshold by setting 
the c and c' parameters to zero in accordance with the expectation that the cross section 
vanishes at the Thomson limit. 

The relativistic Breit-Wigner resonance amplitude Ar[W) for a spin- J resonance R of 
mass rriR is given by 



, 87r(2J+ l)mfl ^V,,,{W)T,,{W)B{r,n^) 
Ar{W) - \ — —2 — — . _-p , (9) 



W ml-W^ -imRTtotiW) 

The energy-dependent total width rtot(W^) is given by 

TtotiW)=Y,^x,xAW) , (10) 

X 

where Xi is tt, K, rj^jj, etc. For J = 2 (the a2(1320) meson), the partial width FxiXal^) 
is parameterized as |23l |: 

( qx{W^) y D2 {qx(W^)rR 
V Qx (mji) J D2 {qx {m%)rR 

where Tr is the total width at the resonance mass, qx{W^) = 
y/{W^ - {mx, + mx.mW^ - {mx, - mx,y)/i2W), D2{x) = 1/(9 ± Sx^ ± x^), 
and TR is an effective interaction radius that varies from 1 GeV"^ to 7 GeV^^ in 
different hadronic reactions [131 ■ For the three-body and the other decay modes, 
r3-body(W^) = FflS(i? — > 3 — body)^^ is used instead of Eq. (fTTI) . All the parameters of 

™R I 1 

the 02(1320) are fixed to the PDG values as listed in Table HVl [21], except for tr, which is 
fitted to be 3.09lEjj^ (GcV/c)-\ consistent with 3.62 ±0.03 (GeV/c)-i determined for the 
72(1270) For the ao(980) and ao(y), the widths are taken to be energy independent. 
For the ao(980), a sirnple Breit-Wigner formula is used instead of the more sophisticated 
formula used in Ref. Ml, i| for the /o(980). This is because the resonance shape appears 
to be symmetric with no indication of the effect of the KK threshold. In fact, a fit with 
the formula in Ref. [l|, Q gives g]^g-/g^T^ — Olg'°^, where qkk {g-KTr) is the coupling of the 
ao(98G) to KK (tttt). 



2. Fitted parameters 



We fit differential cross sections with the parameterized amplitudes for the range 
0.90 GeV < W < 1.46 GeV. There are 19 parameters to be fitted. About one thousand 
sets of randomly generated initial parameters are prepared and fitted using MINUIT [2^ to 
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TABLE V: Fitted parameters. The errors are statistical only. 



Resonance 


Parameter 


Nominal 


M{ao{Y)) fixed 


No ao{Y) 


Unit 


ao(980) 


Mass 
Ttot 

r^^S(r/7i-°) 


982.3lo.? 
75.6 ± 1.6 
128t| 


982.3+0? 
76.911° 
55812 


982.3 ±0.6 
75.6llj 
642 ± 8 


MeV/c^ 
MeV 
eV 




Mass 

Ttot 


1316.8+10 
65.0l^;i 
432 ± 6 


1474.0 (fixed) 

251 + 33 

(11.011^) X 10^ 


(fixed) 


MeV/c^ 
MeV 
eV 



X^/ndf 597.6/429=1.39 704.5/430=1.65 753.6/433=1.74 




0.5 
I cosiJ' I 



0.5 

I COST?' I 



FIG. 11: Differential cross section (da / d\cos6*\ (nb)) (data points) and results of the fit (solid line) 
for the W bins indicated. The dotted, dashed and dot-dashed curves indicate the \S\^ , A-k\DqY2\^ 
and 47r|D2V2^|^ contributions, respectively. 



search for the true minimum and to find any multiple solutions. A unique solution is found 
with /ndf = 597.6/429 = 1.39 {ndf denotes the number of degrees of freedom) for the 
nominal fit, which appears in more than ~ 3% of the cases. The fitted parameters are listed 
in Table IVl The quoted errors are MINOS statistical errors. They are calculated from the 
values obtained by varying each parameter while floating all the other parameters. 



Differential cross sections together with the fltted curves are shown in Fig. [TT] for selected 
W bins. The fit is reasonable as can be seen from these bins and from Fig. [TH where the 
quantities S"^ and D2 are reproduced reasonably well. 

The total cross section (| cos6* \ < 0.8) can be obtained by integrating Eq. ^ as: 

cTtot = O.SS^ + 0A57DI + 0.9831)^ , (12) 
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FIG. 13: Total cross section (|cosS*| < 0.8) and results of the parameterization. Contributions of 

l-DoP and |-D2|^ (not S^, D% and -Di) are also shown. 



where the factors come from the integration of spherical harmonics for | cos0*| < 0.8. The 
measured total cross section is in good agreement with the prediction obtained from the 
sum of the fitted amplitudes as shown in Fig. [I3l 

So far we have used the measured parameter values for the 02 (1320). However, having 
high statistics two-photon production data, one might question their validity, in particular, 
S(a2(1320) — ^ 77). Namely, the determination of this quantity may have been biased in 
past experiments by the presence of non-resonant background, etc. To study this effect, a 
fit is performed where the product V^^B{t]v:^) of the a2(1320) is also floated. The value 
r.y^(a2(1320))S(?77rO) obtained is (145^^^) eV, which corresponds to S(a2(1320) 77) ^ 
(9.4^^;^) X 10-6; it agrees well with the PDG value, (9.4 ± 0.7) x 10"^ within the rather 
large fitting error. Thus, we conclude that the parameters obtained in past measurements 
are reasonable. The large statistical error in our measurement arises because of interference. 

The ao(F) mass is close to that of the a2(1320), which may suggest the possibility that 
the ao(F) is just a contribution of the 02(1320) to the Do wave, which is not taken into 
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account in the above fit. To test this hypothesis, a fit is performed with the ao{Y) removed 
and with Dq and D2 in Eq. ^ replaced by 



where ro2 indicates the fraction of the 02(1320) in the Do wave. We obtain ro2 = l-7lo;4% 
with x^/i^df = 737.7/431. The fit quahty is unacceptably poor. When the ao(l') is restored 
and a fit is performed with Eq. we obtain ro2 = 3Atl;l% with x^/ndf = 580.6/427. 

We conclude that the contribution of the 02(1320) to the Dg wave is small. 

The mass and width of the ao(Y) are significantly smaller than 1474 ± 19 McV/c^ and 
265 ± 13 MeV, the parameters of the ao(1450) in the PDG [2l|. Since the fit without the 
ao{Y) or one where the mass of the ao{Y) fixed to the ao(1450) mass are unacceptable 
(Table [V]), the inclusion of the ao{Y) is required to explain the structure in 5^ near 1.3 GeV 
seen in Fig. [121 Since, the mass and width of the ao(1450) are far from established, we may 
either identify the ao(Y) with the ao(1450) or treat it as another scalar resonance. 



The following sources of systematic errors on the parameters are considered: dependence 
on the fitted region, normalization errors in the differential cross sections, assumptions about 
the background amplitudes, uncertainties from the unfolding procedure, and the uncertain- 
ties of the 02(1320). For each study, a fit is performed allowing all the parameters to fioat 
and the differences of the fitted parameters from the nominal values are quoted as system- 
atic errors. A thousand sets of randomly prepared input parameters are prepared for each 
study and fitted to search for the true minimum and for possible multiple solutions. Unique 
solutions are found very often. 

Two fitting regions are tried: higher (0.94 GeV <W < 1.50 GcV) and lower (0.86 GcV < 
W < 1.42 GeV). Normalization error studies are divided into those from uncertainties of 
overall normalization and those from distortion of the spectra in either | cos 6**1 or W. For 
overall normalization errors, fits are made with two sets of values of differential cross sections 
obtained by multiplying by (1 ± o'e(vK,| cose*|))j where Ce is the relative efficiency error; the 
results are denoted "normalization" errors. For distortion studies, ±4% errors are assigned 
and differential cross sections are distorted by multiplying them by (1 ± 0.1| cos 6**1 ^ 0.04) 
and (l±(VK(GeV)-1.18)/7) (referred to as "bias:| cos 6**1" and "bias:M^" errors, respectively). 

For studies of background (BG) amplitudes, one of the waves is changed to a first- or a 
third-order polynomial. In addition, constant terms are fixed to non-zero values. Uncer- 
tainties from the unfolding procedure are studied by analyzing the differential cross sections 
where a key parameter in unfolding is varied within an allowable range. Finally, the param- 
eters of the 02(1320) are successively varied by their ±la uncertainties. 

The resulting systematic errors are summarized in Table IVTl The total systematic errors 
are calculated by combining the individual errors in quadrature. As can be seen in Table IVH 
the values of T^jB{r]Tr°) for the ao(980) and ao{Y) jump to much larger values in some of 
the systematic variations. This is because destructive interference is preferred in some of 
the studies. A study reveals that for much narrower W regions, e.g. about the ao(980) peak 
region, there exist two solutions corresponding to constructive and destructive interference; 
the latter solution is disfavored in a wider W range in the nominal fit, but favored in some 
of the systematic studies. In this case, taking the sum in quadrature of the individual errors 
will be an overestimation. Thus, we choose the maximum deviation among the different 
systematic variations to estimate the total systematic error. 




(13) 



3. Study of systematic errors 
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TABLE VI: Systematic uncertainties for ao(980) and ao(Y) parameters 





ao(980) 




oourc6 


Mass 


Ttot 




Mass 


Ttot 






(MeV/c^) 


(MeV) 


(eV) 


(MeV/c 


^) (MeV) 


(eV) 


W lit TcingG 


+0.4 
-1.1 


+2.4 
-0.0 


+7.9 
-0.0 


+5.7 
-1.2 


+4.3 
-2.7 


+0.0 
-57.0 


1\ Ul llidllZidLIU-Li 


+0.4 

— 0.6 


+0.2 

—0.2 


+ 17.6 
— 19.5 


+0.9 

— 0.0 


+0.0 
— 1.3 


+20.2 

—46.4 


Riaci" 1 cr\'^ f) 


+0.0 


+0.4 


+0.0 


+3.6 


+6.6 


+0.0 


-0.1 


-0.2 


-9.3 


-0.0 


-18.1 


-179.3 


RiaslV 

1—} 1 tlO , V !■ 


+0.0 
-0.1 


+0.2 
-0.1 


+5.1 
-4.6 


+0.1 
-0.0 


+0.5 
-0.5 


+0.8 
-1.4 


J_) Vjt . J- vc I 1 ^djidjlll. 


+2.9 
— 0.5 


+0.0 
— 1.4 


+0.0 
— 10.6 


+0.0 
— 0.8 


+ 12.2 
— 0.9 


+ 150.2 
— 35.0 


Jj . llii \0 j p dl dill . 


+0.0 
-1.6 


+0.5 
—0.6 


+482.4 
— 6.1 


+17.3 
— 0.0 


+97.8 
— 6.8 


+ 1061.1 
— 141.5 




+0.1 
-0.5 


+0.1 
-0.3 


+0.0 
-9.3 


+0.1 
-0.1 


+0.0 
-0.4 


+0.0 
— 12.0 


BG:Re(Do) param. 


+0.0 
-0.9 


+0.0 
-0.3 


+0.0 
-2.0 


+6.6 
-0.0 


+ 1.9 
-10.2 


+ 18.1 
-215.7 


BG:Im(Do) param. 


+0.0 
-0.9 


+0.0 
-0.3 


+0.0 
-2.0 


+6.6 
-0.0 


+0.0 
-10.3 


+0.0 
-215.7 


BG:Re(7?o) c 7^ 


+0.3 
-0.7 


+0.3 
-0.0 


+0.0 
-9.4 


+6.6 
-0.0 


+0.3 
-14.4 


+0.0 
-229.6 


BG:Re(D2) param. 


+0.0 
-0.1 


+0.5 
-0.0 


+0.0 
-6.3 


+5.4 
-0.0 


+0.0 
-13.4 


+0.0 
-210.8 


BG:Im(L'2) param. 


+0.0 
-0.0 


+0.0 
-0.0 


+0.2 
-3.0 


+0.6 
-0.0 


+0.0 
-1.3 


+0.0 
-18.4 


BG:Re(L>2) c / 


+0.5 
-0.0 


+0.5 
-0.0 


+3.9 
-0.6 


+0.6 
-0.0 


+ 1.2 
-0.3 


+ 11.7 
-6.6 


Unfolding 


+0.7 

-3.9 


+ 17.2 
-9.8 


+501.2 
-31.3 


+9.7 
-3.2 


+0.0 
-7.4 


+0.0 
-213.3 


02: mass 


+0.0 


+0.1 


+0.7 


+0.9 


+0.0 


+0.7 


-0.1 


-0.1 


-0.8 


-0.9 


-0.3 


-5.6 


02: width 


+0.2 
-0.2 


+ 1.0 
-0.8 


+3.9 
-4.1 


+2.7 
-2.8 


+5.4 
-4.7 


+48.0 
-38.9 


02 : r-,^S(r77r'') 


+0.0 

-0.1 


+0.1 
-0.0 


+0.7 
-0.6 


+ 1.1 

-0.6 


+0.7 
-1.6 


+ 12.2 
-25.3 


a2 : rfl 


+0.3 


+0.4 


+5.3 


+0.4 


+2.6 


+8.4 


-0.0 


-0.0 


-0.0 


-0.3 


-1.6 


-10.9 


Total 


+3.1 


+ 17.4 


+501.6 


+24.7 


+99.1 


+ 1073.2 


-4.7 


-10.0 


-43.0 


-4.6 


-32.6 


-255.5 



TABLE VIL Fitted parameters of the ao(980) 



Parameter 


This work 


PDG 


Unit 


Mass 


982.3 


+0.6 +3.1 

-0.7 -4.7 


984.7 ± 1.2 


MeV/c^ 


Ttot 


75.6 ± 1.6 tVnt 


50 - 100 


MeV 


r^-,B(777r'') 


128 


+3 +502 
-2 -43 


240 tfo 


eV 



4- Summary of resonance studies 

To summarize the study presented in this section, once the ampUtudes are parameterized, 
differential cross sections can be fitted to obtain the parameters as described above. Al- 
though the fit itself is not very good as can be seen from /ndf = 1.39, it is stable despite 
the fact that the approach to the minimum is slow; tens of MINUIT runs are needed to 
reach the minimum. The mass, width and T^-yBirji:^) values obtained for the ao(980) and 
ao(F) are summarized and compared to those in the PDG [2l| in Tables [VIll an d lyilll Note 
that the value of the product T^^Bl-qir^) in the PDG [IH is an average of Refs. and (l3l |. 
In both analyses, the total cross section or an event distribution is fitted to an incoherent 
sum of the aq(980) and a2(1320) resonances (with the masses and widths fixed to the earlier 
PDG values [26|) and nonresonant background because of the limited statistics available. In 
our fit, we fully take into account interference among amplitudes. When we follow the same 
procedure as in the previous analyses, which ignored possible interference, we reproduce 
their values with much better statistical errors. 
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TABLE VIII: Fitted parameters of the ao{Y) compared to those of the ao(1450) 



Parameter 


This work 


ao(1450) (PDG) 


Unit 


Mass 


1316.8 in t'f,' 


1474 ± 19 


MeV/c^ 


Ttot 


65.0 til 


265 ± 13 


MeV 




432 ± 6 tllT 


unknown 


eV 



V. ANALYSIS OF THE HIGHER ENERGY REGION 



In this section, we study the angular dependence of the differential cross sections, W de- 
pendence of the total cross section, and the ratio of cross sections for tjtt^ to tt^tt^ production 
in the high energy region, W > 2.4 GeV. 



A. Angular dependence 

As in the analysis of the 7r°7r° process we compare the angular dependence of the 
differential cross sections with the function sm~'^ 9* . A fit with an additional cos^ 6** term 
does not significantly improve the fit quality. Limited statistics prevent us from quantifying 
a possible deviation from the sin~^ 9* behavior when we study the W dependence of the 
data. Here, we only show comparison with a sin~^ 9* parameterization in different W regions 
in Fig. [141 In this figure, the vertical axis is the differential cross section divided by the total 
integral over |cos0*| < 0.8. The curve is 0.602 sin^^0* (not a fit). The numerical factor is 
the differential cross section normalized to a{\ cos9* \ < 0.8). The experimental result shows 
that the agreement is good for W > 2.7 GeV. 



B. Power-law W " dependence 

We fit the W dependence of the total cross section (| cos9*\ < 0.8) in the energy region 
3.1-4.0 GeV, where the lower boundary 3.1 GeV is the same as in the 7r°7r° analysis. The 
fit gives n = 10.5 ± 1.2 ± 0.5, and the corresponding cross section is drawn in Fig. \Wla] 
as well as that of the tt^tt" process in the same angular range. The systematic error is 
obtained from the difference of the central values when we shift the cross section by ±lcr 
at 3.1 GeV and ^la at 4.0 GeV and by factors obtained by connecting Hnearly for W bins 
in between, where a is an energy-dependent part of systematic error at each W point. The 
n value can be compared with n values in other processes that we studied earlier 0, H, Q • 
The results are summarized in Table ITXl from which it is clear that the result for the rjn'^ 
final state is consistent with that for K^Kg (where the fitted W range is wider), but two 
standard deviations higher than that for 7r°7r° . The energy dependence of the latter seems 
to be different from that of the two other purely neutral final states and closer to tt+tt" and 
K~^K~ , although no strict conclusions can be drawn at this level of statistics. 



C. Cross section ratio 



A ratio of cross sections among neutral-pseudoscalar-meson (7r° or 77) pair production in 
two-photon collisions can be predicted relatively easily within a pQCD model. The pQCD 
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TABLE IX: The value of n in atot oc VK " in various reactions fitted in the W and j 008 6**1 ranges 
indicated. 



Process 


n 




W range (GeV) 




1 cos 0* rani 


;e Reference 




10.5 ± 1.2 ±0.5 




3.1 - 4.1 




< 0.8 


This work 


ttV 


8.0 ±0.5 ±0.4 


3.1 


- 4.1 (exclude 3.3 - 


3.6) 


< 0.8 


[9J 


KlKl 


10.5 ±0.6 ±0.5 


2.4 


- 4.0 (exclude 3.3 - 


3.6) 


< 0.6 


15] 




7.9 ± 0.4 ± 1.5 




3.0 - 4.1 




< 0.6 


[3J 


K+R- 


7.3 ± 0.3 ± 1.5 




3.0 - 4.1 




< 0.6 


[3J 



TABLE X: Predictions for the cross section ratio: (T(r;7r")/cr(7r''7r'') in two-photon collisions. 
"Coherent sum" ("Incoherent sum") is the ratio derived from the squared coherent (incoherent) 
sum of the product of the constituent-quark charges. Here, Rf = (/ri//,ro)^, where /,, {fn) is the r] 
{■K°) form factor; the value may be taken to be Rf — 1. The 77 meson is treated as a pure SU(3) 
octet state for the entries in the "octet" row, while "Vp = —18°" is the most probable mixing angle 
between the octet and singlet states from experiments. 



77 in SU(3) 


Coherent sum 


Incoherent sum 


Octet 


0.24/?/ 


0.67i?/ 


Vp = -18° 


0A6Rf 


1.29i?/ 



model in Ref. [lB| predicts the cross section ratio a{ri'K^)/a{TT^TT°) as summarized in Table 1x1 
In the table, Rf = {fril f-^oY , where (/tt) is the 77 (tt^) form factor; the value of Rf is 
not well known and we temporarily assume it to be unity. The ratio of the cross sections 
is proportional to the square of the coherent sum of the product of the quark charges, 
|S]eie2p, in which ei = —62 in the present neutral-meson production cases. We show two 
predictions: a pure flavor SU(3) octet state and a mixture with Vp = —18° for the 77 meson. 
For comparison, we also show in the table the values calculated using an incoherent sum as 
an example of an extreme case. Here, we assume that the quark-quark component of the 
neutral meson wave functions dominates and is much larger than the two-gluon component, 
in obtaining the relations between the cross sections. 



The W dependence of the ratio between the measured cross section integrated over 
I cos 6** I < 0.8 of 77 rjTT^ to 77 7r°7r° is plotted in Fig. fTSlfb). For the 7r°7r° process, 
the contributions from charmonium production are subtracted using a model-dependent as- 
sumption described in Ref. Even though the ratio may have a slight W dependence, we 
average the ratio of the cross sections over the range 3.1 GcV <W< 4.0 GcV as was done 
in other processes for the sake of comparison with QCD and obtain 0.48 ± 0.05 ± 0.04. In 
the averaging, the ratio in the charmonium region (in tt^tt") 3.3 - 3.6 GeV is not used. This 
ratio is in agreement with the QCD prediction if we take Rf — 1. 

D. Comments on charmonium 

It is conjectured that the known cc charmonium states do not decay into the 777r" final 
state with any observable rate, because the 1=1 component should be suppressed. In 
other words, this process could be useful to search for a new charmonium-like particle with 
1 = 1 that would be a candidate for an exotic resonance. In Fig. [161 we show the invariant 
mass distribution of 777r° events with | cos9* \ < 0.4, where a resonance contribution would be 
enhanced. We do not observe any signals of the known Xcj mesons. There is a hint of a peak 
near 3.18 GeV, however, its statistical significance is less than 3a. Therefore, we assumed 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 o 0.2 0.4 0.6 0.8 

|cos e*| |cos e*| |cos e*| 



FIG. 14: The angular dependence of the differential cross sections in different W regions, with the 
normalization to the cross section integrated over |cos0*| < 0.8. The curves are proportional to 
sin~^ 9* and normalized similarly. 




w(GeV) W(GeV) 



FIG. 15: The W dependence of the cross section (|cosS*[ < 0.8) (a). The curve is the power-law 
fit. The cross section of rin^ production is compared to that for n'^n" The W dependence of the 
cross section ratio of rj-K^ to tt^tt" (| cosS*| < 0.8) (b). The line is the average in the 3.1 - 4.0 GeV 
range (the charmonium region, 3.3 - 3.6 GeV, is omitted from the calculation 0|). 

in the above discussion that there are no charmonium contributions in the measurements of 
the cross section described here. 



VI. SUMMARY AND CONCLUSION 

We have measured the process 77 f]n^ using a high-statistics data sample from e+e~ 
colUsions corresponding to an integrated luminosity of 223 fb~^ with the Belle detector at the 
KEKB accelerator. We obtain results for the differential cross sections in the center-of-mass 
energy and polar angle ranges, 0.84 GeV < < 4.0 GeV and | cos^*| < 0.8. 

Differential cross sections are fitted in the energy region 0.90 GeV < W < 1.46 GeV in 
a model where partial waves consist of resonances and smooth backgrounds. The Do wave 
is small, the D2 wave is dominated by the a2(1320) resonance, the S-wave prefers to have 
at least one additional resonance (denoted as ao{Y)) in addition to the ao(980). The mass, 
width and the product r^^B{r]Tr°) for the ao(980) are fitted to be 982.3 tl\ MeV/c^, 



24 



18 

? 16 
S 14 B 

O 

12 

w 

1 1° 

S 8 



4 

E 

5 2 



UUl 








|cos e*|<0.4 



2.8 



3.2 



3.4 
W (GeV) 



3.6 



3.8 



FIG. 16: The W distribution in the charmonium region (| cos^'j < 0.4). 



75.6 ±1.6 



65.0 



+ 17.4 
-10.0 



+3 +502 
-43 



eV and those for the ao{Y) are 1316.8 Z% ifg MeV/c 



MeV, 128 

-I'.l -32.1 MeV, 432 ± 6 ^256^ f'V. The large systematic errors, in particular for two- 
photon widths, originate from an additional solution that favors destructive interference in 
some of systematic studies. The mass and width of the ao(F) are significantly smaller than 
those of the ao(1450). The fact that the obtained aQ{Y) mass is close to the 02(1320) mass 
may suggest that the 02(1320) contribution in the Dq wave is important. However, a fit 
reveals that the fraction of the 02(1320) in the Dq wave is small and it cannot replace the 
ao{Y). Since the fit without it or the fit where its mass is fixed to the ao(1450) mass is 
unacceptable, it is at least a good empirical parameterization. We may still identify the 
ao{Y) with the ao(1450), given that the latter is far from established, or as another new 
scalar meson. We cannot draw a definite conclusion on the existence of the 02(1700). 

The angular distribution of the differential cross sections is close to ~ sin~^ 9* above 
W = 3.1 GeV similarly to the 7r°7r°. In this energy region, the energy dependence of the cross 
section integrated over | cos0*| < 0.8 is well fitted by W~^\ n = 10.5 ± 1.5 ± 0.4, somewhat 
higher (by two standard deviations) than that in the 7r°7r° channel. Although a slight W 
dependence may remain in the ratio, we average the cross section ratio, (T(ry7r°)/cr(7r°7r°) in 
the range 3.1 GcV < < 4.0 GcV and obtain 0.48 ± 0.05 ± 0.04. This ratio is consistent 
with the prediction from a QCD model based on qq production and SU(3) symmetry. 
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